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Our work presents a study on the nonlinear dynamical behavior for a microcavity semiconductor 
containing a quantum well. Using an external periodic perturbation in energy level we observe the 
periodic-doubling, quasiperiodic, and direct route to chaos as forcing strength is changed. For a par- 
ticular case the riddled basin for coexisting periodic and chaotic motions are observed. These results 
suggest that the dynamics of exciton-photon is quite complex in presence of external perturbation. 
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I. INTRODUCTION 



The properties of the light emitted from a quantum well exciton embedded in semi-conductor microcavities have 
drown attention of several theoretical and experimental researchers in the last twenty years [iT-flll] , given their potential 
application in opto-electronic devices [l2l.ll3|. Phenomena which are linked to quantum electrodynamics, such as Rabi 
splitting and Autler-Townes doublets, have been observed [1, EH, Ell- 
in the strong coupling regime when the coupling frequency between the exciton and photon is larger than the 
relaxation frequencies of the medium and the cavity, the vacuum Rabi splitting has been observed. The degeneracy 
between the cavity resonance and the medium is lifted and two lines appear in the intensity spectrum of the exciton- 
cavity system. Semiconductor quantum wells exhibiting narrow absorption lines corresponding to excitonic strengths 
of these resonances make it possible to achieve the strong coupling regime and to observe the vacuum Rabi splitting. In 
the strong coupling regime also, modifications of quantum statistical properties of the emitted light, were theoretical 
predicted and experimentally proved (g, |16H21|. I n 2004, the group of Giacobino proved the existence of the optical 
bistability in the semiconductor microcavities [22j . Bistability occurs when there exist two possible stable states for 
the same set of system parameters. Optical bistability can play a major role for the conception of the optical memories 
paving the way to the realization of the quantum computer. 

Recently, multistability in semiconductor QED systems in the strong coupling regime was observed in different 
experiments [23l [24| . The multistability can arise due to different effect such as polarization-dependent nonlinear 
interactions [23| or thermal induced processes [25[ . 

Either in natural systems or in experimental realizations the presence of an external forcing is unavoidable. Such 
external forcing could be either noise (either from surroundings or inherent within the experiment) or deterministic 
perturbation. The external forcing in nonlinear dynamical systems has been extensively studied. This forcing some- 
times gives important and useful results e.g. stochastic resonances [26|, [27j , chaos control [HI, [2t| , strange nonchaotic 
dynamics [30h32| etc. One aspect for exploring the nonlinear behavior is to scan the parameters space and observe 
how the dynamical complexity depends on the parameters (33M36I ) . For a given set of parameters understanding the 
basins of coexisting attractors are also equally important, particularly, when there are riddled basins [IMi. The 



aim of this Letter is to analyze the dynamical behaviors (particularly routes to chaos [33T - I361 ] and riddling 37-4(|) of 
the intra-cavity photonic intensity and the intensity of the fiorescent light in the presence of periodic signals. 

In this work we explore the dynamics of the field intensities for high excitation regime inside a semiconductor 
microcavity containing a quantum well. We observe that the periodic-doubling, quasiperiodic and direct route chaos 
as forcing strength is being changed. This shows the existence of various types of dynamics over the forcing strengths. 
Furthermore coexisting periodic and chaotic attractors are also observed with riddled basin. 
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The paper is organized as follows. In Sec. II, we review the derivation of exciton-phonon interaction. This is 
followed in Sec. Ill by results and discussion for three particular cases of detuning where various complex dynamics 
are presented. Results are concluded in Sec. IV. 



II. MODEL 



We consider a microcavity, containing a semiconductor quantum well, with width of the order of the wavelength 
A. The quantum well is localized in a position, which corresponds to the maximum of the electromagnetic field. This 
system is described in detail in [4l|. 

Our discussions are limited to the semiconductors with two bands. The electromagnetic field can excite an electron 
of the valence band to the conduction band by creating a hole in the valence band. The electron and the hole interact 
by giving excitonic states, which are similar dependent states to the states of hydrogen atoms. The state IS which is 
the low bound state (fundamental state) has the greatest oscillator strength. It is for this reason that this processing 
of the interaction between exciton - photon often takes into account only this state. Neglecting the effects of the spins 
and the interaction with the phonons we can write a Hamiltonian describing the system |42h44| : 
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The photonic and excitonic modes are quantified according to the normal direction to the microcavity. The vectors 
of waves are thus vectors parallel with the planes of the layers. In a perfect planar cavity, the invariance by translation 
imposes that the exciton with parallel wave vector can couple only with the same photons with parallel wave vector 
k,,=K n . 

The first two terms of the Hamiltonian correspond to the proper energies of photon and exciton, where a k ,/, bx^ 
are respectively the annihilation operators of a photonic and excitonic mode verifying: 
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^>exc,K /r ^ph,k// are the photonic mode and excitonic mode frequencies of the cavity. The third term corresponds 
to the exciton-photon coupling with a coupling constant g. This constant is a function of the oscillator strength by 
unit of area ^, the effective length of the cavity L c , the dielectric permeability e, the mass mo and the free electron 
charge e (SBj : 



9 = V 2rn s f L c S (4) 

and the fourth term represents the saturation of the interaction photon-exciton due to the fermionic nature of the 
electrons and the Holes, with 0, 0, 0, E3| : 
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X x being the radius of Bohr of the bidimensional exciton and it is about 100A in the GaAs microcavity with a 
quantum well. 

The last term describes the interaction between exciton due to the Coulomb interaction between components of an 
exciton and those of a neighboring exciton. The effective potential is given by: 

= Vo" - |^ (6) 

where A is the surface of quantification. We limit ourselves to the case of a pumping having a normal incidence and 
which results in the excitation of only one photonic mode of the cavity with wave vector q = 0. 

The momentum conservation imposes that this photonic mode is coupled only with only one excitonic mode Kit = 0. 
The Hamiltonian describing exciton in the microcavity with a quantum well interacting with the photons is then: 

H = hajp^a^a + haj exc b + b + ihg (a + b — b + aj + — V^b + b + bb 

-lL( a + b + bb + b + b + ab ). (7) 

n sa t 

If the excitation is very weak which correspond to very low density, the distance which separates one exciton from 
its neighbor is very large so that the interaction between excitons is negligible as well as the effects of saturation and 
the exciton can be considered as a pure boson and the exciton-photon interaction can be modeled like an interaction 
between two harmonic oscillators [48145 Oj : 

Hii n = hu! p ha + a + fujj exc b + b + hg [a + b + b + aj . (8) 

The third term of this Hamiltonian indicates the coupling between exciton-photon and represent the exchange 
of excitations. The term a + b describes the creation of a photon and the annihilation of an exciton and the term 
ab + describes the creation of an exciton and the annihilation of a photon. The proper states of this system are called 
polariton states of the cavity. The diagonalization of the Hamiltonian gives energies _E + and E_ of polaritons in the 
cavity: 

E+ = _^LJ_^. ± (9) 

with 



A = ^? + V (10) 
where the difference of frequency between frequencies of exciton and photon: 

<5i = <Jj cav — U! exc . (11) 

The corresponding polariton operators P± are linear combinations of the excitonic and photonic operators: 

P_ = Xb- Ca (12) 

P+=Cb- Xa (13) 
5i + A 



with X and C are the Hopfield coefficients: 
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If we take into account the relaxation phenomena due to the finished lifetime of excitons and photons that can be 
explained by the dissipations of the cavity {"f C av) an d that of the exciton "/ exc . 

We can modalize these effects by introducing the complex energies for the exciton E cav — ihj cav and for the photon 
of the cavity E exc — ihj exc 12]. Then we get complex energies whose real part gives the energy of polariton and the 
imaginary part represents the rates of polariton relaxation: 

T-, E cav -\- E exc ,,lcav lexc , ^A < 1r < 

E± = in ± — — (15) 

2 2 2 y ' 

A 2 = (5 1 -h ( lcav - lexc )f + V. (16) 
The polariton relaxation rates for the upper and lower branches are fH| : 

7+ = Jexc cos 2 (6) + 7 cal) sin 2 (9) 

7- = lexc sin 2 (6) + j cav cos 2 (9) (17) 
where 9 is the coupling angle defined by : 

t 5 (29) = — — Q - - = X. (18) 

£ \0Jcav — OJex) ^1 

In the case where the exciton and the cavity are in the resonance {uj cav = oj e x) dissipations of the two branches are 
equal to the average of excitonic and photonic dissipations: 

7ea:c ^cav 

7+ =7- = 2 ■ ( 9 ) 

If we take into account the effects of nonlinearity due to the excitonic interactions and the effects of saturation we 
must reconsider the fourth and fifth terms of the global Hamiltonian H. 

Let us now derive the evolution of the excitonic and photonic fields in the non linear regime where the density 
of exciton is high. The evolution of the mean field operators in the interaction picture, by considering that the 
fluctuations are weak compared to the average values, can be written as: 



d(a) 
d(b) 

dr ~ 1 " 2 
d(a+) 
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= e(t) + g (b) - k (a) - iA a (a) (20a) 

-9 (a) ~ \ (b) - tA b (b) - 2ia (b+) (b) (b) (20b) 

= e(t) + g (b+) - K (a+) + iA a (a+) (20c) 

(/ : - g (o+) - 1 (6+) + zA b (b+) + 2ia (b+) (b+) (b) (20d) 



dr " v 1 2 

where r is a dimensionless time normalized to the round trip time r c in the cavity: 



_ t 

Tc 

7, k are the dimensionless decay rates of the exciton and cavity photon: 



7 = lex T c ; K = J ph T c 

the nonlinear coupling constant is normalized to i : 
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a = 
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and A a , Aj are the dimensionless detuning 
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As this system is already 4-dimensional and hence the presence of forcing makes the system more complex and 
difficult to understand it analytically. Therefore the numerical simulation is used to uncover the various dynamical 
states using RK4 integrator [5lj . We consider the step size At = 27r/5000 for integration. The dynamical studies are 
studied after removing initial 10 7 data points as transients. We explore here the photon and exciton intensities I a = 
(a + ) (a) and lb = (b + ) {b) inside the cavity. As the fluorescent light is proportional to the mean number of excitons 
(7b), we are also exploring the fluorescent light dynamics. 



There are various physical parameters in this system. We fix the normalized parameters by g = 1.5, K = 0.12, 
7 = 0.015 . These parameters correspond to the experimental values [22[ in units of the inverse of the round-trip 
in the microcavity. We concentrate our analysis to the strong pump field where the non-linear phenomenon exhibits 
interesting dynamics, in this case we choose a normalized amplitude of laser pump e = 200 and the normalized 
excitonic interaction coefficient a — 0.00001. At this set of parameters the system has stable fixed point solution. 
However, particularly e is more prone to external influence. Therefore in this study we consider this parameter only 
to observe the effect of the external forcing. In this work we restrict our analysis to the deterministic periodic forcing 
in e i.e. e : e = 1 + fcos(Qt) where fi is the perturbation frequency while / is the strength of the forcing which is 
considered as bifurcation parameter [52[ . It is clear that for / = the motion is fixed point solution. As we change 
the magnitude of /, various dynamical motions are possible which we will discuss below for different values of the 
detuning A a and A&. We limit our study to the case of the resonance between the exciton and the cavity uj ex = uj p h- 
In the strong regime where the values of the dissipation rates 7 and k are small compared to the coupling constant 
exciton-photon g, the systems exhibits two polaritonic branches at uj ex ±g (see Eq. (15)). We explore the dynamics 
for these two cases where the laser pump is resonant to one of the polaritonic frequencies A a = A& = ±g. We also 
analyze the case of the total resonance between the cavity, laser pump and the exciton A a = A& = 0. 



Let us begin with the case of total resonance A a = A5 = 0. Shown in Fig. HJa) is the spectrum of Lyapunov 
exponents (LE), \ [33j as function of forcing strength /. The variation of LEs indicates the different dynamical 
states as the forcing strength / is changed. Here the zero Lyapunov exponent is plotted as a dotted (blue) line while 
largest two non-zero Lyapunov exponents are plotted as solid (black) and dashed (red) lines. The figures in middle 
row represent the intensity of the photons inside the cavity I a = (a) (black-solid line) and the intracavity excitonic 
intensity = (b + ) (b) (red-dashed line) as function of time at different values of the forcing strengths (b) / = 0.1, (c) 
/ = 0.804 (d) 0.329, and (e) 0.5 which present the different dynamical motions as period-one, period-six, quasiperiodic 
(QP) and chaotic (C) respectively. These different dynamical motions can be further confirmed by Poincare section 
}33j taken at Re < a >= 0, as shown in lowest row. The respective motions are the period-one (single point), period- 
six (six points), quasiperiodic (closed curve) and chaotic (scattered points). These results clearly demonstrate that 
there are various possible dynamical states which depend on the forcing strength. 

As the forcing / increased from zero, the motion goes from periodic (P) to chaotic via QP. The left inset figure 
of Fig. [T] shows the expanded region in near to the transition marked as F\. It clearly indicates that the motion 
goes from periodic to chaotic via quasiperiodic where two Lyapunov exponents are zero. This transition can be 
termed as quasiperiodic route to chaos, a route which has been found is generic in many other dissipative dynamical 
systems |34| - l3a | . Near to the transition F% there is an another route to chaos where the periodic motion gets doubled 
and leads to chaos as the forcing parameter is being increased. This is clearly shown in an expanded Fig. [2] (a) 
which demonstrates that whenever the largest non-zero Lyapunov exponent reaches to the zero value, the period gets 
doubled. This period doubling is further confirmed in Fig. [5] where 100 consecutive maxima of Re < a > (at a fixed 



III. RESULTS AND DISCUSSIONS 



A. 



Case-I: . 



6 



parameter) is plotted as a function of forcing parameter. Here the periods get bifurcated as periods 6 — >• 12 — >• 24.... 
A representative period-six motion is shown in Figs. HJc) and (g). Apart from these transitions there are many other 
periodic windows, e.g. F 4 , F 5 , F7 etc. where one of the above routes occurs. 

However as the forcing strength / is decreased near the marked transition F3, the motion becomes chaotic from 
periodic where there is sudden jump in Lyapunov exponents. This suggests that there is direct transition from periodic 
to chaotic as shown in inset figure near to the transition F3. This suggest that this is an another type of route to 
chaos. Here due to high dimensionality of the system it is difficult to find the normal form to understand the exact 
type of bifurcation. Apart from the above routes/dynamics there are various periodic windows of different periods 
also e.g. near to F4.. These results indicate that in the presence of the external forcing this system exhibits various 
type of dynamics. 

B. Case-II: A a = A b = g 

This situation presents the case where the laser pump is resonant to the upper polaritonic branch. Similar to 
the case-I, we plot in Fig. [3] (a) the largest two Lyapunov exponents, A, as a function of the forcing strength /. 
The figures in middle row show the intensities I a (black-solid line) and lb (red-dashed line) as a function of time 
at different forcing strengths (b) / = 0.06, (c)0.18 and (d) 0.8 for periodic, quasiperiodic, and chaotic attractors 
respectively. These different dynamical motions can be further confirmed by Poincare section taken at Re < a >= 0, 
as shown in lowest row, for respective motions which are the period-one (single point), quasiperiodic (closed curve) 
and chaotic (scattered points) motions. These clearly demonstrate that these various dynamical states also possible 
in this case as well. 

Here as the forcing strength is increased from zero, the period motion becomes quasiperiodic (which is in wider 
region than that of the Case-I) and then is chaotic near to the transition Fx. Here also the transition Fx where 
motion goes from periodic to chaotic is direct in the sense that the periodic motion suddenly becomes chaotic. This 
transition, similar to F3 transition in Case-I, is still not exactly clear due to inability to get its normal form. 

C. Case-Ill: A a = A 6 = -g 

This situation presents the case where we pump resonantly to the lower polaritonic branch. In Fig. 2] (a) are plotted 
the spectrum of Lyapunov exponents, Aj as a function of the forcing strength /. Similar to Figs. [1] and [3] shown in 
middle row are the intensities I a (black-solid line) and lb (red-dashed line) as function of time at different forcing 
strengths (b) / = 0.5 and (d) 1.7728 for periodic and chaotic attractors respectively. The corresponding Poincare 
sections, taken at Re < a >— 0, are shown in lowest row, which confirm the period-one (single point), and chaotic 
(scattered points) motions respectively. These clearly demonstrate that, in this case as well, various dynamical states 
are also possible. 

However the transition from periodic to chaotic and then to periodic is quite different than from the previous 
Cases-I and II. Here we don't find neither period-doubling nor quasiperiodic nor direct route to chaos as observed in 
previous cases. A remarkable feature of this case is that there appears to be a parameter interval, near the transition 
Fx in Fig. Ufa), where the system exhibits wild fluctuations as a function of the forcing strength /. These fluctuations 
persist on small scales, as shown in a blowup of marked box in inset figure. Such fluctuations typically indicate 
co-existing attractors with complicated basins [37T - l40j because of a small change in the parameter /, it can lead to a 
completely different attractor. Here we do observed coexisting periodic (Fig. [SJa)) and chaotic (Fig. HJb)) motions. 
The basin of these multistable attractors are shown Fig. [5jc) with different initial conditions corresponding to periodic 
(open circles) and chaotic (star) motions. In the diagonal region, interwoven initial conditions show the presence of 
complicated basins where both of the attractors coexist. In order to confirm the complicated basins, it is plotted in 
Fig. [5jc) the basin of expanded region of marked box in (d). This demonstrate that basin is riddled |37H4(j | i.e. if an 
initial condition goes to periodic attractor then for any slight change in the / parameter, say within the change <5, it 
may go to chaotic attractor, and vice- versa. However, this phenomenon of riddling can be verified via the computation 
of an uncertainty exponent. Fixing a perturbation, S, and randomly choosing a pair of systems within a region, the 
parameters are termed uncertain if they yield different asymptotic states. The fraction of uncertain parameter pairs, 
q{6), typically decreases with S as a power 

q(S) oc 5 P 

which defines the uncertainty exponent (3 [37j - l40| . Results are plotted in Fig. |SJe) where the exponent is ap- 
proximately zero (the best fit to the data yields f3 — 0.0008 ± 0.0004) which is typical of riddle or riddle-like basins 
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[37H40(. The practical implication is that in this region, the attractor cannot be predicted no matter how small the 
uncertainty in the specification of parameters. This suggests that this system exhibits complicated dynamics which 
not only depends on initial conditions but also on the perturbation parameter. 

IV. CONCLUSION 

In summary, we have studied the dynamical behavior of the photonic and excitonic intensities inside a semiconductor 
microcavity with a quantum well. As the fluorescent light is proportional to the intracavity excitonic intensity, we 
have explored also the fluorescent light dynamics. 

We have observed first that by varying the external perturbation we can change the intracavity field intensities. 
By increasing the amplitude of the external perturbation / and independently of the resonance configurations mainly 
the photon intensity inside the cavity is increased. 

By pumping the upper polaritonic branch (casell) there are more excitons than photons inside the cavity and by 
pumping the lower polaritonic branch (case III) we have the inverse situation (see Figs. 2 and 3). This can be 
interpreted by the fact that excitonic nonlinearity induces more excitons in the upper branch than in the lower one. 
The enhancement of the photon intensity is more pronounced by pumping resonantly to the upper polaritonic branch. 

A modification of the trajectory regularity is observed by increasing the external forcing strength /. As shown in 
Fig 1,2 and 3, by increasing the perturbation amplitude / the dynamics of the system changes from stable point (Fig. 
la, Fig. 2a, Fig. 3a) to chaotic motion (Fig. lc, Fig. 2c and Fig. 3c). The quasiperiodic motion is totally absent in 
the case III. Moreover the basin of attractor in III are quite complicated. 

The Case-Ill shows different behavior in comparison to I and II. The chaos is observed in wide range of parameters 
for I and II. From this we can conclude that the regularity of the dynamics for the polaritonic lower branch is less 
sensitive to the external perturbation. Furthermore we have observed that Hyper chaotic behavior (more that one 
Lyapunov exponents are rather that zero) in I but not in II and III. This means that by pumping in the middle 
between the two polaritonic branches the system presents very complex dynamics. In addition, we can deduce that 
the perturbation affects less the polaritonic lower branch than the upper branch( see Figs. 2 and 3). By pumping with 
a laser frequency detuned from the polaritonic frequencies the system exhibits high sensitivity to the perturbations. 
This effect can originate from the nonlinear effects induced by the exciton-exciton interaction. The nonlinearity affects 
more the frequency ranges detuned from the polaritonic frequencies and induces high dependence of the system to 
the external perturbations. 
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Figure captions: 



Fig.l (Color online) Case-I: A a = Ab = 0. (a) The two largest Lyapunov Exponents with function of forcing 
strength /. Middle row: intensities I a (black-solid line) and lb (red-dashed line) as a function of time at different 
values of the forcing strengths (b) / = 0.1 (period-one motion) (c) / = 0.804 (period-six motion) (d) / = 0.329 
(quasiperiodic motion) and (e) / = 0.5 (chaotic motion). The scale of y-axises in b), c), d), and e) are the same. 
The (f), (g), (h) and (i) are the Poincare section at Re < a >= in the intensities space I a — h for b), c) d) and 
e) respectively. The dotted line in a) represents the zero Lyapunov Exponent. The inset figures are the expanded 
regions near to different transitions as indicated by arrows. 



Fig. 2 (Color online) (a) The two largest Lyapunov Exponents and (b) the bifurcation diagram, Re < a > m (maxima 
of Re < a >) as a function of the forcing strength / near to the marked transition F2 of Fig. 0Ja). The inset figure 
in b) is expanded region of dotted box. 



Fig. 3 (Color online) Case-II: A a = A/, = g. The two largest Lyapunov Exponents as a function of the forcing 
strength /. Middle row: intensities I a (black-solid line) and lb (red-dashed line) as a function of time at different 
forcing strengths (b) / = 0.0.06 (periodic motion) (c) / = 0.18 (quasiperiodic motion) and (d) / = 0.8 (chaotic 
motion). The (e), (f) and (g) are the Poincare section at Re < a >= in intensities space I a — lb for b), c) and d) 
respectively. The dotted line in a) represents the zero Lyapunov Exponent. 



Fig. 4 (Color online) Case-Ill: A Q = Af, = —g. The largest two Lyapunov Exponents as a function of the forcing 
strength /. Middle row: I a (black-solid line) and h (red-dashed line) as a function of time at forcing strengths 
(b) / = 0.5 (periodic motion) and (c) / = 1.7728 (chaotic motion). The (d) and (e) are the Poincare section at 
Re < a >= in intensity space I a — lb for b) and c) respectively. The dashed line in a) corresponds to the zero 
Lyapunov Exponents. The inset figures in a) and e) are the expanded regions of marked boxes in the respective figures 
with more data points. 



Fig 5 (Color online) Case-Ill: A Q = A& = —g. The presence of coexisting attractors (a) period and (b) chaotic 
as function of time at a forcing strength / = 1.77273. (c) Then basin of these coexisting attractors in Re < a > 
vs Im < a > with fixed Re < b >= Im < b >= 0. (d) is an expended basins of marked box in c). (e) shows the 
fraction of uncertain parameter pairs (out of a sample of 300) as a function of the parameter perturbation S (see text 
for detail). 
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